Population of closed-channel molecules in trapped Fermi gases with broad Feshbach resonances 
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We compute the fraction of closed-channel molecules in trapped atomic Fermi gases, over the entire range 
of accessible fields and temperatures. We use a two-channel model of BCS-Bose-Einstein condensation (BEC) 
crossover theory at general temperature T, and show that this fraction provides a measure of the T dependent 
pairing gap. Our calculations, containing no free parameters, are in good quantitative agreement with recent 
low T measurements in 6 Li. We present readily testable predictions for the dependencies of the closed-channel 
fraction on temperature and Fermi momentum. 



PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 

The superfluid phases of ultracold, trapped atomic Fermi 
gases present an exciting opportunity to study how superfluid- 
ity is changed as the system evolves from the weak attraction 
limit of Bardeen-Cooper-Schrieffer (BCS) state to the strong 
attraction Bose-Einstein condensation (BEC) regime 1 1J. Fes- 
hbach resonances arise when a bound state ("closed channel 
molecules") of a two-body spin-singlet potential lies near the 
continuum onset of a spin-triplet scattering state ("open chan- 
nel"). As a consequence of an applied magnetic field B, prox- 
imity to these resonances can profoundly change the two-body 
(s-wave) scattering lengths a. The well studied Feshbach res- 
onances in both 6 Li and 40 K are very wide as a function of 
B, compared to the Fermi energy. Moreover, the focus is on 
unitary scattering where a diverges. For these reasons one of- 
ten resorts to a one-channel version of BCS-BEC crossover 
theory although attention has also been paid to the two- 
channel variant f° r wide as well as narrow Fes- 
hbach resonances. 

In recent experiments 01 on the broad resonance in 6 Li 
at low temperatures, T, a laser drives transitions between 
dressed molecules (hybridized with open channel atom pairs) 
and an excited molecular singlet, in the process removing ex- 
cited atoms from the trap. The exponentially decreasing num- 
ber of remaining trapped atoms as a function of probe du- 
ration time is used to determine that fraction Z of dressed 
molecules/pairs corresponding to the closed channel. Here 
we present the analogous theoretical calculations of Z, as a 
function of a and for general T. We demonstrate very good 
agreement with experiment, and provide predictions for future 
finite T experiments. The small size of this fraction means 
that for many properties the two-channel model reduces to the 
single-channel description. We emphasize however, that the 
two-channel model contains more information making it ca- 
pable of addressing a wider range of physical systems and 
experiments, and of testing the fundamentals of BCS-BEC 
crossover theory. 

Our theoretical approach reduces to the generalized 
Leggett-BCS theory [ 1 ] at T = 0. At finite T, the noncon- 
densed as well as condensed contributions enter in the ap- 
propriate combinations, so that a future measurement of the 
total closed-channel fraction will provide a measure of the In- 
dependent pairing gap A(T). This gap is, in general, distinct 
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from the superfluid order parameter, A sc (T) , except at T = 0. 
Finally, we present results for the closed-channel fraction as a 
function of the global Fermi wavevector hp, at different fields 
B. At unitarity, the dependence is a simple proportionality. 
When two-channel physics is important, as it must be in the 
determination of the closed-channel fraction, then the energy 
scale which characterizes the width of the Feshbach resonance 
also enters into the problem. In this way this fraction is non- 
universal and assumes different values for 6 Li and 40 K. 

We begin with a discussion of the superfluid phase. A gas 
of Fermi atoms in the presence of a Feshbach resonance con- 
tains open-channel fermionic atoms as well as closed-channel 
molecules. Because of coupling g between the open and 
closed channels the three propagators for the open-channel 
fermion pairs [t pg (Q)], the closed-channel molecules [D(Q)] 
and the single (open-channel) fermion states [G(K)] are all 
highly interconnected. Here and throughout, we use a four- 
momentum notation: K = (ia> n ,k) and Q = (iSl m ,q), 
J2k = fe^EnEk' an d tne i r analytical continuation, 
ioj n — > ui + i0 + ,ifl m — > ft + i0 + , where u> n — (2n + 
VjirhsT /h, f2 n = 2mrkBT/h are the odd and even Matsub- 
ara frequencies. The T-matrix scheme we employ to treat 
these coupled propagators is derived from the equations of 
motion for the Green's functions, and, importantly, it natu- 
rally leads to the self-consistency conditions of the standard 
T = mean field theory JlJ. It can be shown ||J that 
the pairs are described by the pair susceptibility x(Q) = 
Eif G o(Q ~ K)G(K)tpl_, 2 where G de P ends on a B CS- 
like self energy E(isT) » -A 2 G a (-K )<p£. Here G (K) is 
the noninteracting fermion Green's function. Throughout this 
paper ip^ = exp{— k 2 /2k 2 ,} introduces a momentum cutoff, 
where fcrj, represents the inverse range of interaction, which is 
assumed infinite for a contact interaction. 

Our Hamiltonian is the standard boson-fermion model 
1I3H5 HI0I1 in which there are only fermion-boson and fermion- 
fermion interactions. Central to our analysis at finite T are 
noncondensed pairs J3l which we characterize by a T-matrix 
t Pg (Q) = U eff (Q)/[l + U eff (Q)x{Q)], where U eff is the 
effective pairing interaction which involves the direct two- 
body interaction U as well as virtual excitation processes as- 
sociated with the Feshbach resonance At small Q, t pg 
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can be expanded as 



tpg{Q) 
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The parameters appearing in Eq.^are discussed in more de- 
tail in Ref. 11(11 . Here Z^ 1 is a residue and Q, q the pair 
dispersion. The latter parameter as well as the effective pair 
chemical potential /J, pa i r both depend on the important, but 
unknown, gap parameter A through the fermion self energy 
S. The decay width Tq is negligibly small for small Q below 
T c . We caution here that the inverse residue "Z" appearing 
in Eq. Q is not the same as the quantity Z used in Ref. |8]. 
More technical details about the various residues can be found 
in Ref. full . 

To determine A, one imposes the BEC-like constraint 
li pair = which yields tz}(Q -> 0) = = U^JQ) + X (0), 



i.e., 



rr-i rn\ , 1 ~ 2 /(- g k) ,„2 



(2) 



so that A formally satisfies the usual BCS gap equation 
with quasiparticle dispersion — (ik — /i) 2 + A 2 ip 2 ,, 
where et = H 2 k 2 /2m is the fermion kinetic energy, /i is the 
fermionic chemical potential, and f(x) is the Fermi distribu- 
tion function. 

Physically, one should view A as reflecting the presence of 
bosonic degrees of freedom. In the fermionic regime (/i > 0), 
it represents the energy required to break the pairs. It can 
be seen |9] that A contains contributions from both noncon- 
densed and condensed pairs, whose densities are proportional 
to A 2 g (T) and A 2 C (T), respectively. One has a constraint on 
the number of pairs H 1 Oil which can be viewed as analogous to 
the usual BEC number constraint 



A 2 (T) = A 2 C (T) + A 2 9 (T) 



(3) 



The total contribution of noncondensed pairs is readily com- 
puted in terms of A via 
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In analogy with the standard derivation of BEC, one can then 
compute the number of condensed pairs associated with A sc . 

Fermions are the fundamental particles in this system, so 
that their chemical potential [i is determined from the number 
conservation constraint 



ri/ + 2n b0 + 2n b = n f + 2n\ ot . 



(5) 



where n&o and rib represent the density of condensed and 
noncondensed closed-channel molecules, respectively, n£°* is 
their sum, and n/ = 2 J^k ^(-^0 * s tne atom i c density asso- 
ciated with the open-channel fermions. Here 



n b = -2_j D (Q) ~ Zb X! fo (^q - Mpoir) 

Q q 



where b(x) is the Bose distribution function. The renormal- 
ized closed-channel molecule propagator D(Q) is given by 
the same equation as Eq. Q with a different residue Z^ 1 — > 
Z b . We compute A (and /i) via Eqs. and l|5}, to determine 
the contribution from the condensate A sc via Eqs. (0 and @. 
Note that at T = 0, A pg = so that all pairs are condensed as 
is consistent with the mean-field BCS-Leggett ground state. 

At T — 0, a non-zero excitation gap for all k constrains n / 
so that one can identify our closed-channel fraction 2n t b ot /n — 
2nbo /n with the quantity Z in experiment. In this picture the 
ground state consists of n/2 pairs with fraction 2rn,v/n in the 
closed channel. A short range attractive interaction (appro- 
priate for cold Fermi gases), leads to a picture in which all 
fermions are paired, since the attraction extends to the entire 
Fermi sphere. Importantly, this is consistent, at least at unitar- 
ity, with the observed exponential behavior of the remaining 
total particle number in the trap found in Ref. |8j Thus, we 
presume here that the ground state has no unpaired fermions. 

The interaction parameters U and g which appear in the 
Hamiltonian can be related to their experimental counterparts 
Uq and go. The latter are, in turn, determined by the open- 
channel background scattering length a^g and the Feshbach 
resonance width W: Uq = Anabgh 2 /m = 8wkFabg(EF / kp) 
and <7q = |C/oW|- Here Ep = h 2 kp/2m is the global nonin- 
teracting Fermi energy. From the Lippmann-Schwinger equa- 
tion Jill we obtain l/t/ = 1/U — l/U c , where 1/U C = 
~ Ek(^k/^ e k) ' s m e value of the interaction correspond- 
ing to unitary scattering (a = oo). Similarly, we define 
U* = 4ira*H 2 /m and the "experimental" magnetic detuning 
2/q such that 



U* = Un 
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As in Ref-HH defining T = 1/(1+U /U C ), one has U = TU , 
g = Tgo and v = vq — Tg^/U c . In order to make contact with 
experiment, it is convenient to define a two-body counterpart 
of U* in Eq. Q: 



U 2 b = 4TraH 2 /m = U - gl/vo ■ 



(8) 



which diverges at vq = 0. Hereto = (B— Bo)Afi a , where Bq 
is the resonance field, and AfjP — 2^lb for 6 Li is the differ- 
ence in the magnetic moment between open-channel pairs and 
the closed-channel molecules 1 12J. Here /is is the Bohr mag- 
neton. For the wide Feshbach resonances in 40 K and 6 Li, the 
dimensionless parameters 1/kpa and 1/kpa* are very close 
to each other for the magnetic fields addressed in Ref. 1 8]. 

It is useful, now to rewrite the closed-channel density in 
terms of more physically accessible parameters. Building on 
Refs. j^lS' one can show 01 ti&o oc A 2 c . It follows from our 
T ^ formalism 1 10] that when a condensate is present 
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where we have used some simple algebra to rewrite Z g 



(6) 9 /[(2a* — v)XJ + g ] . Importantly, in this way we can con- 
clude from Eq. that the total closed-channel contribution 
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1: (Color online) Calculated closed-channel fraction 
/N = 2N b0 /N at T = as a function of magnetic field 
B (main figure) at Tf = 0.2/iK for 6 Li in a harmonic trap. The 
inset shows 2Nto/N as a function of l/fc^a (black curve) at Tf = 
0.4/iK. The red circles are experimental measurements from Ref. 1 8], 
and the blue diamonds are theoretically calculated data using spe- 
cific values of Tf from experiment 11311 . which vary between 0.18- 
0.662 fjK. 



i%l ot — Z g A 2 . The simplicity of this last result reflects the 
fact that D(Q) and t pg (Q) [of Eq. Q] share the same denom- 
inator. The relative probabilities for the (hybridized) pairs to 
live in the closed and open channels are essentially fixed for 
all T < T c and given by Z b = Z g /Z and 1 — Z b , respectively. 

To represent the trap, we use the local density approxima- 
tion (LDA) by replacing // — > fi(r) = \i — V(r). Here /i is 
the global chemical potential and V(r) = m,ui 2 r 2 /2 for a har- 
monic trap with angular frequency ui. We solve Eqs. @ and 
(|2} at each r for given /i and then self-consistently adjust /x to 
satisfy the total number constraint N = J d 3 rn(r). We de- 
fine Nbo = j d 3 rn b o(r) and similarly for the noncondensed 
molecules N b , and the total iV b * ot = N b0 + N b . Here N w 
represents the trap average of the superfluid order parameter. 
Moreover, this trap average is proportional to the value of the 
order parameter at the center so that N b o/N cx A^fr = 0). 
Thus, this full two-channel calculation provides a theoreti- 
cal underpinning for the simple interpretation, provided in 
Ref. 1 8], of their experiments. 

In the unitary regime where \a bg /a*\ <C 1, Eq. @ implies 
ribo ~ ^sc/.9o' at eacn P om t r. Therefore, we can deduce that 
at unitarity 
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At unitarity, 2N b0 /N scales with hp ■ Here we have used the 

relationship N = n(0) J d 3 r [n(r)/n(0)] and n(0) cx E 3 J 2 . 
The fraction in front of \[Ef / ' g\ is dimensionless and inde- 
pendent of Tp. 

In Fig. [2 we plot the calculated molecular fraction 2N b o /N 
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FIG. 2: (Color online) Closed-channel fraction as a function of T/T c 
at unitarity for Tf = 0.4/iK for 6 Li in a harmonic trap. The black, 
red, and blue curves are the condensed (2Nw>/N), noncondensed 
(2N b /N) and total (2AT* ot /jV) fractions, respectively. Here T c = 
0.273T F . In the BEC regime 2Nl ot /N is relatively T independent 
below T c . 



at T = for 6 Li (black curve and blue diamonds) as a func- 
tion of the magnetic field B as compared with experimental 
data (red circles) from Ref. |8]. The experimental values for 
Tp vary from one data point to another. Our calculations at 
each field value use as input the specific experimental value 
for Tp, plotted as blue diamonds, in the figure and inset. How- 
ever for the continuous curve plotted in the main figure, we 
use Tp = 0.2/iK which better represents the values in the 
BCS regime II 311 . In the inset we used Tp = 0.4/iK, since 
this best reflects the average value over the entire range of 
data points. These slightly different choices for Tp reflect 
the fact that the two plots have different horizontal axes, and, 
therefore, amplify small errors in different ways. As will be 
shown below, in the BCS regime the results are particularly 
sensitive to Tp. We take the parameters for the Feshbach res- 
onance from Ref. (3: W = 300 G and a field-dependent 
a bg = o° ff /[l + <*(B - B )}, where a° bg = -1405a , 
a = 0.0004 G^ 1 , B = 834.15 G and a = 0.529A is 
the Bohr radius. For Tp = 0.4/iK, this uniquely determines 
U = -5.90E F /k F and g = -771E F /k 3 J 2 at unitarity. 

In Fig.^the roughly a factor of 3 difference between the- 
ory and experiment in the BCS regime may in part be asso- 
ciated with the observed deviation from an exponential time 
dependence |8]. Due to the extremely small gap (< T) at 
these high fields, unpaired fermions may be thermally excited. 
However, it should be noted that this particular thermal "cor- 
rection" alone would be in the wrong direction (see below). 
This may point to the need for a more systematic treatment of 
dynamical effects in this regime. Overall, agreement with the 
data is quantitatively very good 1 8] everywhere but in the BCS 
regime; there are no adjustable parameters flUl . 

We turn to the T dependence of the closed-channel fraction, 
which is plotted at unitarity in Fig. [2] This fraction was shown 
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FIG. 3: (Color online) Predicted closed-channel fraction 2Nbo/N at 
T = as a function of \JTf oc k,F at different fields for 6 Li in a 
harmonic trap. At unitarity (B — 834G, black line), 2Nto /N oc kp. 
On the BCS side (B = 900G, red, and B = 1000G, blue), this 
becomes a higher power law. On the BEC side (B < 834G, see 
the inset, B increases from top to bottom), the fraction becomes less 
sensitive to Tf, and is determined by two-body physics in the deep 
BEC limit. 



more readily the closed-channel molecules decay into the 
open channel so that the steady state molecular fraction re- 
mains small, as observed experimentally. We have shown here 
that measurements of the closed-channel fraction will com- 
plement other techniques for obtaining the pairing gap. While 
broad Feshbach resonances often lead to universal behavior at 
unitarity 11711 . a single dimensionless parameter kpa is inad- 
equate for determining quantities such as 2Nj,o/N, which is 
intrinsically associated with two-channel physics. 

We thank R. Hulet for sharing data, and J.E. Thomas, C.A. 
Regal, M. Greiner, and C. Chin for discussions and the OCTS 
workshop. This work was supported by NSF-MRSEC Grant 
No. DMR-02 13745 and by the Institute for Theoretical Sci- 
ences and DOE, No. W-31-109-ENG-38 (QC). 

Note added. — After submission of this Letter, two papers 
appeared on the T = closed-channel fraction within a ho- 
mogeneous system based on a different model for the ground 
state II 811 and on a limiting case 1 19] of the present theory 
in which direct interfermion interactions and closed-channel 
self-energy) are ignored. 



above to be proportional to the pairing gap squared A 2 (T). 
The condensed fraction 2 N^o /N (black curve) decreases as T 
increases from zero and vanishes at T c . At the same time, the 
noncondensed fraction 2Nb/N (red) increases from zero, as 
determined by Eq. 0. It has a maximum slightly below T c . 
The total fraction 2Nl ot /N (blue) decreases monotonically 
with T. Note that it decreases very slowly at low T. This 
justifies our comparison in Fig.^between the T — calcula- 
tions and the low T measurements. Using Eqs. (|9jl and i ll Oi l, 
and in conjunction with the temperature insensitivity of the 
various residues (Z g and Z^), one can measure the T depen- 
dence of the pairing gap, presuming that T is determined from 
isentropic sweeps 1 16]. The sensitivity of the total fraction to 
temperature increases with field. 

Next we explore the relationship between the closed- 
channel fraction and the Fermi temperature Tp. In Fig. [5] we 
plot 2N b0 /N at T = as a function of y/Tp ozkp for differ- 
ent magnetic fields, from the weak pairing BCS (B = 1000G, 
blue curve) to the strong pairing BEC (B = 700G, cyan 
curve). It can be seen that the plots reflect our earlier theo- 
retical prediction that 2Nbo/N oc kp at unitarity. As the field 
increases, 2Nbo/N varies as a higher power of kp. In con- 
trast, as the field decreases in the BEC regime, it becomes less 
sensitive to kp. For 6 Li, it approaches unity in the deep BEC 
limit, where it is dominated by two-body physics. More gen- 
erally, 2Nf > o/N increases faster (slower) on the BCS (BEC) 
side than the simple proportionality found at unitarity, as is 
confirmed by Fig. [5] 

In summary, we find that the wider is the resonance the 
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